United States Patent im 

Wu et al. 


[ii] Patent Number: 
[45] Date of Patent: 


[54] PIPELINED DIGITAL SAR AZIMUTH 
CORRELATOR USING HYBRID 
FFT/TRANSVERSAL FILTER 

[75] Inventors: Chialin Wu, Pasadena; Kuang Y. Liu, 
Los Angeles, both of Calif. 

[73] Assignee: The United States of America as 

represented by the Administrator of 
the National Aeronautics and Space 
Administration, Washington, D.C. 

[21] Appl. No.: 314,928 

[22] Filed: Oct. 26, 1981 

[51] Int. C1.3 GO IS 9/00 

[52] U.S. Q 343/5 CM; 343/5 FT; 

343/5 DP 

[58] Field of Search 343/5 CM, 5 DP, 5 FT, 

343/9 PS, 5 CF; 364/728, 726, 724 

[56] References Cited 

U.S. PATENT DOCUMENTS 


4,003,054 1/1977 Goldstone 343/5 FT X 

4,034,370 7/1977 Mims 343/5 DP X 

4,045,795 8/1977 Fletcher et al 343/5 DP X 

4,132,989 1/1979 Arens 364/728 X 

4,292,634 9/1981 Wu et al 343/5 CM 

4,298,950 11/1981 Smith, Jr 364/726 

4,355,311 10/1982 Jam 343/9 PS X 


OTHER PUBLICATIONS 
Lawrence R. Rabiner & Bernard Gold, Theory and 


Application of Digital Signal Processing Prentice Hall 
1975 Chapter 6 pp. 356-437. 

Primary Examiner — Maynard R. Wilbur 
Assistant Examiner — Brian Steinberger 
Attorney, Agent, or Firm — Paul F. McCaul; Thomas H. 
Jones; John R. Manning 

[57] ABSTRACT 

A synthetic aperture radar system (SAR) having a 
range correlator (10) is provided with a hybrid azimuth 
correlator (12) for correlation utilizing a block-pipe- 
lined Fast Fourier Transform (I2a) having a predeter- 
mined FFT transform size with delay elements (Z) for 
so delaying SAR range correlated data as to embed in 
the Fourier transform operation a corner-turning func- 
tion as the range correlated SAR data is converted from 
the time domain to a frequency domain. A transversal 
filter (126) connected to receive the SAR data in the 
frequency domain, and from a generator (146) a range 
migration compensation function, D, to a programma- 
ble shift register (30) for accurate range migration com- 
pensation; weights, W ;, to multipliers (35-38) for inter- 
polation, and an azimuth reference function, <f>j, in the 
frequency domain to a multiplier 42 for correlation of 
the SAR data. Following the transversal filter is a 
block-pipelined inverse FFT (12c) used to restore azi- 
muth correlated data in the frequency domain to the 
time domain for imaging. The FFT transform size is 
selected to accommodate the different SAR azimuth 
aperture lengths, number of looks and prefiltering re- 
quirements. 

14 Claims, 11 Drawing Figures 
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PIPELINED DIGITAL SAR AZIMUTH 
CORRELATOR USING HYBRID 
FFT/TRANSVERSAL FILTER 

5 

ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 1Q 
Aeronautics and Space Act of 1958, Public Law 85-568 
(72 Stat. 435; 42 USC 2457). 

BACKGROUND OF THE INVENTION 

This invention relates to digital signal processing for 15 
synthetic aperture radar (SAR), and more particularly 
to pipelined digital signal processing for producing 
real-time high-resolution SAR images using a hybrid 
processor comprised of a fast Fourier transform (FFT) 
for correlation in the SAR azimuth (along-track) direc- 20 
tion and a time-domain transversal filter to accommo- 
date the dispersion of the SAR response in the range 
direction. This dispersion is caused by the variation of 
target distance to the SAR transmitter-receiver over the 
SAR aperture length, and is often referred to as the 25 
range migration effect that must somehow be compen- 
sated in the SAR data processing. 

The utility of synthetic aperture radar, under investi- 
gation for approximately three decades, has been found 
to be a useful instrument for a number of geological and 30 
oceanographical applications. The launch of a space- 
borne SAR on the SEASAT satallite in June, 1978, 
marked the advance of SAR remote sensing technology 
from a conventional airborne environment into an earth 
orbiting spacebome environment. 35 

Two important aspects of SAR systems are the high 
data acquisition rate and the complicated data reduction 
process. Real-time data reduction to form images for 
the airborne SAR system was accomplished in the mid 
1970’s. Real-time spacebome SAR processing for pro- 
ducing high quality images is still a continuing chal- 
lenge. Much of the operational SEASAT SAR image 
formation process is currently being handled by an 
optical processor and a moderate throughput software- 45 
based digital SAR processing system. In order to meet 
the objectives of potential candidates of future opera- 
tionally oriented spacebome SAR missions in which the 
image data will be used in a prompt fashion to derive 
surface state of temporal change data, a real-time or 50 
near real-time SAR processing capability is required. 

Techniques for achieving this real-time processing to 
form SAR images involve a number of engineering 
disciplines. A survey of the state-of-the-art SAR pro- 
cessing indicates the following conclusions: (1) the the- 55 
oretical framework for SAR sensor operation and data 
processing is reasonably well understood; (2) many 
techniques for SAR data reduction for forming images 
have been devised, but their performance characteris- 
tics and design tradeoffs have not been fully established; go 
and (3) the major cost factor of a high quality real-time 
SAR processor is the need for a tremendous number of 
high speed electronic devices and associated intercon- 
nective networks. 

Existing candidate architectures for a real-time SAR 65 
processor include the time-domain transversal filter and 
the frequency-domain multistage Fast Fourier Trans- 
form (FFT). A comparison of these two types of ap- 
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proaches yields differences in the arithmetic computa- 
tion and processing control requirements. 

In digital processing of SAR data, the raw data that 
has been quadrature demodulated and digitized is first 
range compressed. This range compression is required 
only if the radar employs a coded waveform in its trans- 
mitted pulses. The output of this range correlation pro- 
cess is sequential in the range dimension. Therefore, 
before starting a two-dimensional azimuth correlation, 
one needs to perform the so-called corner-turning or 
matrix transpose on the range correlated data. This is a 
quite time-consuming process. If two-dimensional azi- 
muth correlation is done entirely in the time domain, 
then fine interpolation between data at discrete loca- 
tions must be performed to trace the maximum ampli- 
tude response. Also the azimuth dimension is usually 
quite large, so a large time-domain transversal filter is 
required to perform the azimuth correlation. All of 
these considerations make the time-domain approach 
undesirable. 

But the alternative to perform a two-dimensional 
azimuth correlation entirely in the frequency domain is 
also undesirable. For severe range migration, this alter- 
native requires a large amount of hardware to do the 
two-dimensional transform. Furthermore, this alterna- 
tive requires the transform of the two-dimensional ref- 
erence function. Hence changing the reference function 
is not easy. 

What is required is a digital synthetic aperture radar 
(SAR) azimuth correlator which produces real-time, 
high-resolution SAR imagery utilizing the frequency 
domain for range and azimuth correlation without re- 
quiring excessively large memory capacity for the SAR 
data being processed, and a comer-turn technique that 
will alleviate the complicated matrix transpose problem 
in a two-dimensional azimuth correlation following the 
range correlation with spectral filtering performed in 
the azimuth dimension, and a range transversal filtering 
in the time domain to conveniently accommodate the 
range migration effect. These features yield on efficient 
and modular SAR processor with which an exact SAR 
correlation algorithm may be implemented. 

SUMMARY OF THE INVENTION 

In accordance with the present invention, an SAR 
azimuth correlator uses a combination of frequency 
domain and time domain processing for the two-dimen- 
sional azimuth correlation. The azimuth correlator is 
comprised of a block-pipelined FFT, a transversal filter, 
an azimuth reference multiplier, a block-pipelined in- 
verse FFT and a generator for range migration com- 
pensation and azimuth reference functions used in the 
transversal filter. The range correlated data, received 
by the block-pipelined FFT sequentially line by line, is 
delayed properly in each stage of a block-pipelined 
FFT. Thus corner-turning operation is embedded in the 
FFT processor. Target response interpolations and 
range migration corrections in the range dimension are 
performed in an interpolation filter. This interpolation 
filter is followed by an azimuth reference multiplier. 
Since the interpolation weights are fixed and the azi- 
muth reference function only change gradually with 
range, the control part of the time-domain processing is 
very simple. In case of severe range migration, namely 
dispersion of the SAR response in the range direction 
caused by the variation of target distance to the SAR 
transmitter-receiver over the SAR aperture length, 
often referred to as the range migration effect, the inter- 
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polation weights and azimuth reference will be com- 
bined to perform the true hybrid azimuth correlation 
algorithm (frequency domain in azimuth dimension and 
time domain in the transversal filter for the range migra- 
tion compensation). The output of the reference multi- 5 
plier is fed into the block-pipelined inverse FIT. The 
magnitude of the inverse FFT output is then generated 
(squared) for imaging. Thus SAR images are formed 
sequentially line-by-line in a continuous manner at the 
output of the SAR processor. The magnitude of single- 10 
look images will also be used to provide clutterlock and 
autofocusing adjustments for the azimuth portion of the 
control processor. 

The novel features that are considered characteristic 
of this invention are set forth with particularity in the 15 
appended claims. The invention will best be understood 
from the following description when read in connection 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 2Q 

FIG. 1 is a system block diagram of a hybrid (time- 
domain and frequency-domain) SAR processor accord- 
ing to the present invention. 

FIG. 2 is a block diagram of a pipelined FFT for 
range correlation. 25 

FIG. 3 is a block diagram of a pipelined inverse FFT 
for range correlation. 

FIG. 4 is a switch timing diagram for the pipelined 
FFT of FIG. 2, and the inverse pipelined FFT shown in 
FIG. 3. 30 

FIG. 5 is a timing diagram for the hybrid azimuth 
correlator of FIG. 1. 

FIG. 6 is a functional block diagram of a block-pipe- 
lined FFT for the azimuth correlator of FIG. 1. 

FIG. 7 is a functional block diagram of an exemplary 35 
transversal filter used to implement the range migration 
compensation and range interpolations and azimuth 
reference multiplication function in the system of FIG. 

1 . 

FIG. 8 is a functional block diagram of a block-pipe- 40 
lined inverse FFT for the azimuth correlator of FIG. 1. 

FIG. 9 is a diagram useful in understanding the cor- 
ner-turn operation of the azimuth correlator. 

FIG. 10 is a flow diagram of an 8-point block-pipe- 
lined FFT or an di33 8 data block for corner turning. 45 

FIG. 11 is a block diagram illustrating the organiza- 
tion of a range migration compensation function and 
azimuth reference function generator in the system 
shown in FIG. 1. 

DESCRIPTION OF PREFERRED 50 

EMBODIMENTS 

The SAR signal processing system disclosed in FIG. 

1 is based on a hybrid SAR azimuth correlation algo- 
rithm developed recently to handle SEAS AT SAR 55 
type of signal processing problems which are character- 
ized by a large correlation compression ratio and a 
severe range curvature. It is comprised of a range corre- 
lator 10 and an azimuth correlator 12, together with a 
programmed digital control processor 14 which pro- 60 
vides in the usual manner functions indicated by the 
functional blocks 14a, b and c labeled RANGE REF- 
ERENCE FUNCTION GENERATOR, RANGE 
MIGRATION COMPENSATION FUNCTION 
AND AZIMUTH REFERENCE FUNCTION GEN- 65 
ERATOR and TIMING CONTROL CIRCUIT. The 
range correlator, which compresses the transmitted 
pulse waveform into a narrow pulse, provides fast fre- 
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quency-domain correlation in the range dimension 
using a pipelined Fast Fourier Transform (FFT) 10a, a 
digital multiplier 106 and a pipelined inverse Fast Fou- 
rier Transform (FFT) 10c. The azimuth correlator per- 
forms the synthetic aperture processing with a combina- 
tion of fast frequency-domain correlation in the azimuth 
dimension and a time-domain convolution in the range 
dimension to accommodate the two-dimensional range 
curvature effect. The control processor provides the 
range and azimuth reference functions and timing con- 
trol to the rest of the SAR processor. 

The significance of the design is a block-pipelined 
Fast Fourier Transform (FFT) 12a in the azimuth cor- 
relator which incorporates both the transform process- 
ing and the memory comer-turn (matrix-transpose) 
operations. Subsequently, a transversal filter 126 pro- 
vides range migration compensation in the azimuth 
spectral domain. A block-pipelined inverse Fast Fourier 
Transform (FFT) 12c restores the SAR data to the time 
domain for imaging through a squarer 12c?, and option- 
ally a multiple-look processing circuit as disclosed in 
U.S. Pat. No. 4,292,634 granted to the Administrator of 
the National Aeronautics and Space Administration 
before this application was filed by the Administration. 
Spacecraft data carrying the system, namely position, 
velocity and acceleration vectors, are inputs to the con- 
trol processor which utilizes the predetermined range 
pulse waveform for generating the range reference 
function and the spacecraft data (position, velocity and 
acceleration vectors) for generating the range migration 
compensation and the azimuth reference functions. In- 
terpolation weights generated by the azimuth portion of 
a control processor are set inputs to multipliers of the 
transversal filter in a conventional manner as will be 
described more fully hereinafter. The detected single 
look range data are also fed into the azimuth portion of 
the control processor to provide a closed-loop clutter- 
lock and autofocusing adjustment. 

Before proceeding with a more detailed description 
of the invention shown in FIG. 1, a brief review of the 
hybrid correlation algorithm will first be given. The 
major emphasis of this disclosure will then be to focus 
on the new SAR correlation architecture. 

Review of the Hybrid Correlation Algorithm 

In a typical SAR operation, phase encoded pulses are 
transmitted from the radar antenna in a direction per- 
pendicular to the radar flight path, and echo signals are 
received in time-sequenced bursts. It is conventional to 
represent the received echo signals in a two-dimen- 
sional form with time delay within an echo burst period 
as the vertical axis of the presentation and along-track 
elapsed time as the other axis. After designating the two 
spatial variables r and x for the cross-track (range) and 
along-track (azimuth) distances to the radar sensor, it 
may be shown that the echo response s(x,r) of an ex- 
tended target which has a radar cross section cr(x,r) can 
be written as: 

s(x, r) — cr(x, r) (g)h(x, r) ( 1 ) 

where ® stands for a two-dimensional convolution, and 
h(x,r) is the impulse response of the SAR and is ex- 
pressed as: 

h(x,r)=h\ (x, r)®hi(x, r) 


( 2 ) 
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The hi and I 12 functions are the azimuth and range re- 
sponses, respectively. Equations 1 and 2 show that a 
reconstruction of the target reflectivity o-(x,r) can be 
performed by a correlation of the SAR echo data s(x,r) 
with range response h 2 (x,r) and subsequently with a 5 
modified azimuth response, which is a convolution of 
the hi function and the compressed range waveform. 
Filtering for range response h 2 (x,r) is a standard pulse 
compression process in radar system and can be per- 
formed by several different types of pulse compression 10 
schemes. The emphasis of SAR processing is upon azi- 
muth correlation to compress the dispersed azimuth 
waveform into a narrow one. 

If linear FM (frequency modulation) phase encoding 
is used in the transmitted waveform, then range com- 15 
pression will produce a “sine x over x” sin x/x type of 
radar response in the range dimension. Correlation of 
this range sin x/x response and the azimuth response 
hi(x,r) is given by the following equation: 


20 


h\(x, r) = WJ.x) exp 


. 4tt n(x) "I 

-J 1 J «('- 


( 3 ) 


nW) 


where 

/= radar wavelength 

W a (x)=SAR antenna response in azimuth 
ri(x)= range delay function between sensor and tar- 
get 

8(x)= Dirac’s Delta function 
results in a waveform: 


25 


30 


h{x, r, /■„) = X g,{x, r„) S(r - di) 
1=1 


( 4 ) 35 


where 


giix r a ) = W a (x) 


sin C 0 (dj - n(x, '•<;)) 
Cddi - n(x, r 0 )) 


exp-AK*. r„) 


( 5 ) 


40 


<M*> r 0 ) = 


4v r t (x, r 0 ) 


( 6 ) 


Equation 7, and performing Fourier transforms on the 
resulting equation with respect to x yields: 


0(«. r 0 ) = 2 S(u, di + r 0 ) Gf(u, r 0 ) 
1= 1 


( 8 ) 


45 


C 0 is a factor depending on the range resolution, and d, 
is the range of the i-th range bin. Note that a parameter, 
r 0 , which denotes the target slant range at the center of 
the SAR antenna illumination, is used to reflect the 
range dependent nature of the range delay function r/(x) 50 
in Equation 3. Also note: that in Equation 4, the range- 
correlated azimuth response function is expressed as a 
summation of azimuth response along slant range bins. 
The extent of the summation indicates the amount of 
range migration which traverses through n range sam- 55 
pie elements. Each g,(x,r 0 ) is a one-dimensional wave- 
form. 

The gi functions for different i have overlapped am- 
plitude response due to the sin x/x function. However, 
a continuous phase function given by <f(x,r 0 ) * s pro- 60 
vided. 

The reconstructed reflectivity <r(x,r 0 ) is obtained 
from: 


<Hx,r 0 )=f J i(a+x, 0 + r o ) hi f (a,0,r o ) dad (3 


( 7 ) 65 


where s(x,r 0 ) is the range-correlated SAR echo, and 
means complex conjugate. Substituting Equation 4 into 


where fl(u,r 0 ), S(u,d,-f r 0 ), and G,*(u,r 0 ) represents the 
Fourier transforms of cr(x,r 0 ), s(x,d,-j-r 0 ), and g/*(cr,r 0 ), 
respectively. Equation 8 implies that a filtered image 
line for targets at range distance r 0 is a linear superposi- 
tion of several filtered azimuth line spectra at range 
distance from r<,+d/to r 0 +d«. This filtering approach is 
indeed a combination of frequency-domain processing 
in azimuth and a time-domain convolution which in- 
volves reference multiplication and coherent summa- 
tion in the range dimension. Also, note that this correla- 
tion approach as defined in Equation 8 is exact for a 
target response which is assumed to be invariant in the 
azimuth dimension. The filter function G/* can be up- 
dated periodically along the range dimension to accom- 
modate the range dependent nature of the Doppler 
frequency rate. The approach described in Equation 8 is 
referred to as the hybrid correlation algorithm for SAR 
processing. 

This hybrid approach is applicable for any general 
two-dimensional correlation problem, where the corre- 
lation reference is invariant over the two-dimensional 
data block. For SAR processing with small to moderate 
range curvature or migration effect, the dispersion of 
SAR response is small. This hybrid approach has ad- 
vantages over two-dimensional FFT fast correlation in 
two ways: (1) reference function can be continuously 
updated as processing proceeds in range swath direc- 
tion; and (2) a short convolutional filter in range dimen- 
sion may have a higher computational efficiency than 
the FFT fast correlation implementation. A special 
processor which is based on this hybrid algorithm and 
lends itself to a pipelined operation will now be de- 
scribed. 

Review of One-Dimensional Pipelined FFT and Its 
Application to SAR Range Correlation 

The range correlator 10 is a straightforward match 
filter comprised of a pipelined FFT processor 10a 
shown in FIG. 2 followed by a pair of range reference 
function multipliers 10 b, one for each of two outputs C 
and D of the FFT, and a pipelined inverse FFT proces- 
sor 10c shown in FIG. 3 which receives inputs E and F 
from the multiplier 108. The range reference function is 
the complex conjugate of the transmitted baseband 
signal. 

Referring now to FIG. 2, an N-point decimation-in- 
frequency (DIF) pipelined FFT is shown for range 
correlation, where N equals 2" and n is an integer. There 
are n FFT butterfly modules consisting of a delay ele- 
ment 20 at one of two inputs of a butterfly logic struc- 
ture 22, and a delay element 24 connected to the one of 
two outputs of the butterfly logic structure, to which a 
“twiddle factor” is applied from a table TFT in a read- 
only memory. This twiddle factor is a weighting value 
W ; which may be precomputed and stored for the FFT 
process in a conventional manner. The delay elements 
maybe shift registers implemented with first-in-first-out 
(FIFO) memories of N/2 stages for the delay element 
20, and N/4 stages for the delay element 24. The nota- 
tion Z-W2 and Z~ N/4 thus means a delay of N/2 and a 
delay of N/4 for the two delay elements 20^and 24, 
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respectively. Other FFT butterfly modules are similarly by the authors. What is unique in the present application 

implemented with the pair of delay elements having half of a pipelined FFT is that the input signal is applied to 

the s ages of the corresponding pair of the preceding both inputs of the first butterfly module, to one input 

module, except the last which has only one delay ele- through a delay and to the other input directly, 

ment of one stage. 5 ’ . , „ , ... TT - ™ , 

For correlating long sequences, the overlap-save Azlmuth Correlator Arch itectu re Using Block 

correlation scheme is used. This scheme requires that Pipelined PM 

the input of the pipelined FFT processor 10a (a digi- A key element of the invention is the block-pipelined 
tized and quadrature demodulated SAR echo) be de- FFT processor 12a which accepts sequential radar echo 

layed before entering one input of the first butterfly unit 10 samples and performs a Fourier transform in the azi- 
22 and that the same input data be fed directly into muth dimension. The output of the block FFT proces- 

another input of the butterfly unit 22. The data is then sor consists of sequential data points corresponding to a 

processed through n stages for N data points. Since the column of Fourier coefficients in an azimuth trans- 

overlap-and-save scheme produces N/2 bad data points, formed two-dimensional data block, 

only one output on the pipelined inverse FFT 10c 15 The data handling procedure is described in the fol- 
shown in FIG. 3 is used as the input to the azimuth lowing formulation. Let {s,y> be a two-dimensional 

correlator. The data are preferably complex logarith- (range and azimuth) block of SAR echo data consisting 

mic numbers to reduce the size of the delay elements of N range bins (di points) and M echo lines (d 2 points), 

and butterfly logic structure. The corresponding azimuth transformed data are de- 

To implement this overlap-save correlation scheme 20 noted by {S,y}. The sequence of sampled range echo 
each module is connected to the next by a switch that data input to the block FFT processor is according to {. 

alternately connects one module to the next directly and . . SijSj+ lj , . . snj,sij+ i, . . . } for N samples per echo line 

then crosses the connections, with suitable delays in one (column) of data. The’output data points of the block- 

input and one crossed output of each module. There are pipelined FFT processor are also in a sequence of in- 

therefore n switches (SW) for n FFT butterfly modules. 25 crea sing i. Let {G*,y} be the filter coefficients and {%> 
The timing for the first four of these switches shown in the filtered spectra. For each column of FFT output 

FIG. 4 can easily be generated by a sampling clock the hybrid correlation algorithm calculates {SliJ} 

counter having stages 1 through n numbered from the of Equation 8 according to: 
most significant stage having a value 2" to the least 
significant stage having a value 2 l . The outputs of the 30 

stages thus alternate the switches at different time inter- n i7 -= 2 Si + k~\,j- G‘ k . ' ' 

vals, with SW2 operating at twice the rate of SW1, SW3 k= 1 ’ J 

at twice the rate of SW2, and so forth, so that as the first 

module stores N/2 points and then switches, the next Equation 9 above indicates a n-pomt convolver which 
stores N/4 points and then switches, and so forth. 35 °P erates on the Founer coefficients inside each column 
The pipelined inverse FFT processor 10c shown in sequence and then outputs a filtered spectra {fly}. 
FIG. 3 is organized in much the same way, but re- This time-domain convolution accomplishes the con- 
versed, with two inputs E, F and one output G. Note relation function for SAR response that is dispersed in 

that the delay elements of delay times Z -1 through range dimension. A block-pipelined inverse FFT pro- 

Z-W4 are in an order inversed from that of the forward 40 cessor 12c inputs the filtered spectra {fly} and outputs 
FFT processor shown in FIG. 2, and that for timing the SAR correlated echo signals {Ay} in range sample 
switches are numbered inversely from n-1 through 1. sequence. 

Consequently the timing waveforms shown in FIG. 4 As noted hereinbefore, the range correlated data is 
are applied to the switches in the inverse order, as indi- fed into the azimuth correlator sequentially, line-by-line 
cated by the notation in FIG. 4, with the switch nearest 45 along the range dimension, and it should be recalled 
the input having the shortest switching period and the that the hybrid azimuth correlation algorithm requires 
switch next to the output having the longest switching performing a Fourier transform on the results of range 

cycle. Also note that the first butterfly module has no correlation along the azimuth direction. Therefore a 

input delay element, and that as the delay elements get corner-turn operation is needed to line up data points 
progressively shorter in the pipelined FFT processor 50 along the azimuth direction. 

10a, they get progressively longer in the inverse pipe- A novel technique is to implement the corner-turn 
lined FFT module 10c. operations among the FFT stages. To illustrate this 

Although, for convenience, only four stages (butter- technique, consider the d t Xd 2 data block shown in 
fly modules) are shown for the pipelined FFT proces- FIG. 9, where di and d 2 are the number of data samples 

sor, it should be understood that the number of stages is 55 along the range direction and the azimuth direction, 
in fact n for an N-point decimation-in-frequency (DIF) respectively. A d 2 point block-pipelined decimation-in- 
pipelined FFT and inverse FFT. For a further discus- frequency (DIF) FFT 12a for a diXd 2 data block is 

sion on the organization and operation of such a for- shown in FIG. 6. For correlating these long sequences, 

ward and inverse FFT, reference is made to Theory and the overlap-save correlation scheme is used. As in the 

Application of Digital Signal Processing by Lawrence R. 60 pipelined FFT for correlation in the range dimension, 
Rabiner and Bernard Gold, published by Prentice Hall, this scheme requires that the input data be delayed 

Inc. (1975), and in particular chapter 6 titled “Spectrum before entering one input of the first butterfly unit and 

Analysis and the Fast Fourier Transform” and chapter the same input data is fed directly into the other input of 

10 titled “Special-Purpose Hardware for the FFT” the first butterfly unit. Thus the input data are over- 

wherein, among other things, the authors introduce and 65 lapped by half a block. 

use the now accepted term “twiddle factors.” FIGS. 2 Note that the delay on the first FFT stage is equal to 
and 3 for the forward and inverse pipelined FFT are half the data block long. This is because in the first stage 

functional block diagrams using the symbols employed of the FFT, the two data points which should be oper- 
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ated on along the azimuth direction are separated by 
half a block. In the second stage of the FFT, the input 
is delayed by i block, i.e., delay time Z~ rf l rf 2 /4 , to line 
up the data points. The rest of the FFT stages are oper- 
ated similarly. The flow diagram of an 8-point block- 5 
pipelined FFT for a di X 8 data block is shown in FIG. 

10. Note that each point in FIG. 10 represents a line of 
di samples along the range direction, and not a data 
point in an 8-point DIF pipelined FFT. This is impor- 
tant to the present invention for it is here that the cor- 10 
ner-tum operation takes place in the azimuth correlator. 

The procedure thus outlined is shown with a timing 
diagram in FIG. 5. Note that range correlated SAR 
data in range sample sequence is provided as continuous 
input data to the forward block-pipelined FFT where 15 
transformation of overlapped blocks takes place. A pair 
of n-point convolution filters, each comprised of a trans- 
versal filter which applies range interpolation weights 
to the data, as shown in FIG. 7, and an azimuth refer- 
ence multiplier, also as shown in FIG. 7, the inverse 20 
block-pipelined FFT receives the outputs of these con- 
volution filters and processes it still in overlapped 
blocks. From the inverse FFT, only one output is used, 
namely the output L shown in FIG. 8. As a conse- 
quence, only alternate ones of the two overlapped 25 
blocks of data are used, as though a switch, SW shown 
in FIG. 5, is alternated between output terminals of the 
inverse FFT. These blocks are merged into a continu- 
ous stream which, upon squaring, provides azimuth 
correlated SAR image data in proper range pixel se- 30 
quence. 

The design of the block-pipelined FFT processor 12a 
is analogous to a conventional one-dimensional pipe- 
lined FFT processor as shown in FIGS. 2 and 3. Note 
that the block-pipelined FFT processor 12a for the 35 
azimuth correlator performs transformation only along 
the azimuth dimension. By substituting each point stor- 
age or delay in the conventional one-dimensional FFT 
processor with a shift-register storage of one line (or 
column) of data which consists of N sampled data 40 
points, the transform manipulation is performed on data 
separated by N-l points in the input sequence. These 
points correspond to a row of data along the azimuth 
dimension in a two-dimensional presentation of SAR 
data shown in FIG. 9. The output of the block-pipelined 45 
FFT processor now represents a column of spectral 
coefficients denoted by {Sy}. The idea of extending one 
point to one line shift delay in the block-pipelined FFT 
processing is analogous to the extension of a one- 
dimensional discrete convolver into a convolver with 50 
tapped registers for a time-domain SAR azimuth corre- 
lator. 

It should be noted that the Fourier coefficients after 
the transformation in the azimuth dimension are output 
in range column sequence, i.e., {. . Sij,Si+ij, . . }. Hence, 55 
the filter multiplication as expressed by Equation 9 can 
be performed by the transversal filter 12 6 implemented 
with a n-point one-dimensional transversal filter. The 
phases of Gkj are approximately constant for a given j 
which refers to a constant azimuth position for most of 60 
the SAR operation of the transversal filter 126 in nomi- 
nal situations. A simplification is to separate G/ Ci/ into the 
amplitude weights (W), phase (<j>j) and a coarse delay 
for range migration compensation introduced through a 
programmable shift register as shown in FIG. 7, which 65 
illustrates the organization of the transversal filter 126 
of FIG. 1. The reference coefficient {G*y} can be grad- 
ually updated in range by the control processor to ac- 
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commodate the range dependent nature of the target 
response. The coefficients are changed for every col- 
umn of spectral coefficients. 

The inverse transform of the spectrally filtered data is 
performed by the block-pipelined inverse FFT proces- 
sor 12c, which is essentially another block FFT proces- 
sor having its architecture inverse to that of the forward 
block-pipelined FFT processor 12a. With bit-reversed 
sequence of the transformed coefficients as input, the 
sequence of the inversely transformed output will be 
back to the original echo line order in time-domain. 

Using the FFT fast-correlation approach as incorpo- 
rated in the hybrid algorithm per each block of transfor- 
mation, the number of valid output points equals the 
difference between the sizes of the transform block and 
that of the reference function. It is for this reason that 
the transformation must be processed with successive 
data blocks overlapped by at least: the size of the corre- 
lation reference function as shown in FIG. 5. In order to 
complete the processing in one pass, the block-pipelined 
FFT processor 12a must produce transformed blocks 
that are based on overlapped data in time-domain. 

The block-pipelined inverse FFT processor 12c is 
operated in inversed manner of the forward transform 
processor 12a. It has dual inputs corresponding to the 
dual outputs of the forward FFT. The two output ports 
of the last stage of the block-pipelined inverse FFT 
processor 12c shown in FIG. 8 are data in the time 
domain. Only one output, delayed by one-half of the 
transform block size with respect to the other output, is 
used. A connection is. made at the input of the magni- 
tude generator (squarer 12d), as shown in FIG. 5, to 
accept one-half of each data block at the output of the 
final two ports to produce valid azimuth correlator 
output in a continuous fashion. Of course, the correla- 
tion reference dimension must be less than half of the 
transform block size. 

Note that since the correlator both accepts and out- 
puts data in the form of sequential range samples, the 
internal block processing operation for the transform 
operation in the correlator is entirely transparent to the 
imaging operation. This is extremely desirable to facili- 
tate the line-scan recording to produce a continuous 
long strip of image. The magnitude generator following 
the inverse transform processor is provided to produce 
image data for visual display. 

From the foregoing, it is evident that the azimuth 
correlator is comprised of a block-pipelined FFT fol- 
lowed by two transversal filters for range migration 
corrections, two azimuth reference multipliers and a 
block-pipelined inverse FFT. Two transversal filters 
and two azimuth reference multipliers are required 
because the block-pipelined FFT processor has two 
outputs, output i (H) and output 2 (1) as shown in FIG. 6, 
each of which must be processed separately through a 
transversal filter. Consequently, the structure of the 
transversal filter shown in FIG. 7 for the output H is 
repeated for the output I. Both outputs J and K of these 
two structures are inputs to the block-pipelined inverse 
FFT as shown in FIG. 8. The magnitude of the inverse 
FFT output is then computed to produce SAR imaging 
data. 

Major issues in assessing the flexibility of a SAR 
correlator include the change of the reference function 
width (the synthetic aperture width) and number of 
looks. As mentioned above, the size of the range and 
azimuth reference function cannot be greater than one- 
half of the transformation block due to the 50% over- 
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lapping in the transform processor. When a processor is 
designed to handle the largest possible synthetic aper- 
ture width, it naturally accommodates any shorter syn- 
thetic aperture length without requiring reprogram- 
ming of the correlator. The throughput is either in real- 5 
time rate or in the rate allowed by the limit of the inter- 
nal clock speed of the processor. 

Using the radix-2 pipelined FFT design, the proces- 
sor can accommodate 2 l looks of SAR processing where 
1 is a non-negative integer. The multiple look processing *0 
is accomplished by a subdivision of the azimuth spectra. 
Each look is produced from a segment corresponding to 
a 2 _ / part of the azimuth spectrum. After the azimuth 
spectra are multiplied by the corresponding filter coeffi- 
cients, the products are fed into the inverse block FFT 15 
processor, bypassing the first 1 stages of the radiax-2 
FFT operators consisting of both of the butterfly unit 
and the delay memory for that stage. It should be kept 
in mind that as a characteristic of the FFT operation, 
adjacent Fourier coefficients of a 2~ ! part of the spec- 20 
tral are output in a sequence separated by 21 — 1 data 
intervals. With the first 1 stages bypassed, the inverse 
transform size becomes 2 m ~ I , where 2 m =M is the origi- 
nal transform block. The switching period for the data ^ 
selector following the last stage must be reduced ac- 
cordingly to maintain a continuous flow of valid output 
points. Also, overlay delay memory is needed to pro- 
duce a multiple look image. 

Modular expansion of the FFT processor to accom- 3Q 
modate a longer synthetic aperture length is straight- 
forward due to the inherent modularity of the pipelined 
FFT processor architecture. 

The operation of the pipelined processor can be syn- 
chronized by one clock timing signal source. This facili- 35 
tates the use as a real-time SAR processor, because the 
processor clock frequency can be a multiple of the radar 
pulse-repetition-frequency which often varies accord- 
ing to the sensor flight conditions. 

SAR processing to produce 25 m resolution, 4-look 40 
SEAS AT SAR images requires a synthetic aperture 
length of approximately 1024 radar transmitted pulses 
or echo lines per look. At an optimized sampling fre- 
quency, the number of range samples over the 100 km 
swath width is approximately four thousand. The trans- 45 
form block size of the azimuth correlator is set to be 
2048. The forward block FFT processor requires a 
memory of approximately 4 K X 1.5x2048= 12 Mega 
words, and the inverse block FFT processor requires 4 
Kx2 K=8 Mega words. Assuming 32 bits are used in 50 
each comples work the total memory required is ap- 
proximately 20 Mega, 32-bit words. The inverse block 
FFT processor performs, a four, 512-element transfor- 
mation for the four look processing using the scheme 
described in the previous section. A data throughout 55 
rate of 10 MHz is required to produce images at a real- 
time rate. 

The vast amount of memory required represents one 
major driving factor for the cost of implementing such 
a design. Another major bottleneck of implementation 60 
is the lack of the state-of-the art serial memory of large 
capacity for azimuth correlation. To remedy this, one 
can use the parallel implementation of multiple block 
FFT processors. Each processor accommodates a part 
of the swath width. This subswath partitioning may also 65 
be required to reduce the average data rate so that con- 
ventionally available random access memories (RAM’s) 
and multipliers can readily be used. 
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Two transversal filters are needed, as noted hereinbe- 
fore, one for each output H and I of the block-pipelined 
FFT. The block diagram of the transversal filter is 
shown in FIG. 7. To correct for range migration, a 
coarse delay is provided by a programmable shift regis- 
ter 30 under control of the control processor at the 
input of the transversal filter for range migration com- 
pensation. This coarse delay provides the alignment of 
the first target response along each line. Subsequent 
data samples along each line are delayed by the same 
amount. The coarse delay is changed on a line-to-line 
basis. The target response of the same range distance 
will automatically be lined up along the azimuth direc- 
tion by the coarse delay. 

It should be noted that the timing of the digital sam- 
pling pulses is such that they form a rectangular grid in 
the two-dimensional range and azimuth plane. How- 
ever, the range delay history of a point target is a two- 
dimensional curve as shown in FIG. 9. Therefore, the 
discrete digital radar echo samples do not coincide with 
the loci of the maximum return from a point target. The 
use of a nearest-neighbor approximation on the curved 
target delay history will cause a loss in the pixel signal- 
to-noise ratio as well as a broadening of the target re- 
sponse to occur. Hence range interpolation to introduce 
points nearer to the target return curves is desired for 
achieving a better performance. 

The range interpolation is achieved by a weighted 
sum of four range migration corrected samples as 
shown in FIG. 7, i.e., the output x(t) of the adder 40 is 
given by: 

4 (10) 

*(') = .2, WtxfyiT) 

i=l 

where X,(nT) are the outputs of a tapped delay line 
comprised of shift registers 31-34 and W, are fixed inter- 
polation weights of the form: 


Sin^-frJ (r — iT) 

(•f-) (/ - iT) 

In the above equations T is the digital sampling time, 
and the multiplications are carried out in multipliers 
35-38. The sum is then formed in an adder 40. The 
output of the interpolation filter is resampled and fed 
into an azimuth reference multiplier 42. The purpose of 
the azimuth reference multiplier shown in FIG. 7 is to 
provide necessary phase compensation. 

The block diagram of a d2-point DIF block-pipelined 
inverse FFT with corner-turn memory for a di X d2 data 
block, where d2=2 m , is shown in FIG. 6. The inverse 
FFT is similar to the forward FFT described above, 
except here different delays are used among the FFT 
stages, and only one output of the last FFT stage is used 
due to the fact that the overlap-and-save scheme pro- 
duces \ block of bad data. Thus, the output of the in- 
verse FFT is sequential line-by-line along the range 
direction. This output is fed into the squarer 12 d (FIG. 
1) which provides an output to an imaging device and to 
the range migration and azimuth reference function 
generator 146 for autofocusing and clutterlock adjust- 
ments. 



4 , 471,357 


13 

The organization of the generator 14 b (FIG. 1) for 
the range migration compensation and azimuth refer- 
ence functions will now be described with reference to 
FIG. 11. While reference will be made to discrete func- 
tional blocks which may be implemented with hard- 5 
ware, it should be understood that in practice the func- 
tions indicated by the blocks are carried out by a pro- 
grammed control processor. 

The generator 146 accepts the estimates of Doppler 
center frequency f o and frequency rate fb from an out- io 
side source, either a control microcomputer or an oper- 
ator keyboard, and generates a set of initial azimuth 
references. The initial references are used in azimuth 
correlation, and the output of the azimuth correlator is 
fed back to the azimuth reference function generator for 15 
adjustments in the (d and fp estimates. A new set of 
reference functions is generated, and the feedback con- 
tinues. The generator is therefore comprised of an azi- 
muth reference phase generator section 44 for generat- 
ing range migration and azimuth reference functions 20 
(coefficients), and an interpolation weights and range 
migration correction delay generator section 45, both 
sections operating on the best estimates of fp and fp, 
which are continually revised by a clutterlock section 
46 and an autofocusing section 48. The clutterlock sec- 
tion provides an error signal Afp to refine the estimate 
of Doppler frequency fp based on the single-look image 
data from the azimuth correlator 12 (FIG. 1). The au- 
tofpcusing section provides an error signal Afp to refine 
the estimate of Doppler frequency shift rate fp based on 39 
the single-look image data. The functional operations of 
these three sections will now be described in order. 

Section 44 generates the azimuth phase <f) over the 
S AR aperture duration. That phase follows a quadratic 
function of the elapsed time along track, t, and is given 35 
by the following equation: 

<S(t)=fbt+ Jt'z>t 2 (12) 
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input is integrated by the accumulator 44 b and the sum- 
ming circuit 44a. The sum of (5 and integrated a is 
further integrated by a summing circuit 44 d and accu- 
mulator 44c to produce the value of <J>(u). This value is 
used in a lookup table to determine the complex azi- 
muth reference coefficients applied to the multipliers in 
the transversal filter 126 of the azimuth correlator as 
shown in FIG. 7. 

A similar operation is carried out by elements 45a, b, 
c, d and e to obtain a value <f>(t) along track with inputs 
a' and {S' computed in the arithmetic unit from its inputs 
as follows: 

=h(fD+Ci^fo) (12) 

P =/c2(/b+C24/b) (18) 

where k2 is a constant determined by the radar fre- 
quency pulse repetition rate. The phase <f>(t) is divided 
by a first constant bi in a divider 45/ to obtain a value of 
range migration on the order of meters. This range 
migration value is used in a lookup table 45g to deter- 
mine the interpolation weights Wi, W2, W3 and W4 
applied to the transversal filter as shown in FIG. 7. The 
25 range migration value out of the divider 45/ is also 
applied to another divider 456 where it is divided by 
another constant b2 to produce a delay value D applied 
to the programmable shift register 30 in the transversal 
filter shown in FIG. 7. The constants bi and b2 are 
scaling constants for the range interpolation and sample 
positions determined from the radar transmitter and 
range sampling frequencies, which are generally related 
by an integer multiplier. 

The clutterlock section 46 computes the mean image 
energy Ei, E2, E3 and E4 of four looks over the same 
target area and then obtains the error signal, Afo, ac- 
cording to the equation: 


A Fourier transform of the complex phase function 
e >‘ HO results in the following form: 

(i3) 

H(u ) = /VKO) = Ae 1 2/D 


r 

where A is a constant. The increment of u in equally 
spaced frequency domain is a constant. Therefore, an 
array of the phase values of H(u) for a fixed increment 
in u can be generated by a double integrator comprised 
of elements 44a, b, c, d and e. More explicitly, the phase jq 
of H(u) is 


4>(u) = fa 


, Id , 1 

+ 7“ “ + -rr 

jd Aid 
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The a and fi inputs to the summing circuits 44a to 44c 
are quantities derived by an arithmetic unit 50 from its 
inputs as follows: 


a — k\ 


1 


Id + CiAfo 


(15) 60 


Id + C-jLfp 
2 (Id + QA/d) 


(16) 
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where ki is a constant determined by the radar pulse 
repetition frequency, Ci and C2 are constants adjustable 
to obtain the best feedback loop performance. The a 


iF D =(El+Ez-E3-E4) (19) 

Each of the values E are produced by integrating 
successive single looks. The E/ values all correspond to 
the same target area. Therefore the measurement used 
as the error signal Afo is independent of any target 
scene variation. An architecture for computing this 
measurement is disclosed in FIG. 8 of a copending 
application Ser. No. 969,761 filed by Chialin Wu and 
Vance C. Tyree on Dec. 12, 1978. This error signal Afo 
is integrated and added to the Doppler frequency esti- 
mate fo by the arithmetic unit. 

The autofocusing section 48 computes the crosscor- 
relation of the first and last look images which corre- 
spond to the leading and trailing edge of the synthetic 
aperture. Consequently, a running succession of looks 
are stored in a data storage unit 48a for correlation of 
each one entered with the corresponding last or trailing 
edge look. Assuming the ith pixel in azimuth of the first 
look is denoted as P//, and the corresponding pixel of the 
last (4th) look to be overlaid to P// is P4; (assuming 4 
looks total), the crosscorrelation derived error signal is 
given by the following equation 


A/b =(/*l.iX A./- 1 +P\,iX Pj.i-2) 

-(fi,;Xi > 4,;+i+t’i,iXP4,i + 2) (20) 

where the bar over the products refers to the mean of 
the quantity below. This error signal Afo is integrated 
by a summing circuit 48c and accumulator 48 d. 
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It should be noted that the error signal Afois updated 
periodically through a switch 48a, i.e., the error signal is 
com’ uted for a period T and then applied to the sum- 
ming circuit and accumulator at the end of each period 
by effectively closing the switch momentarily. The 5 
integrated error signal from the accumulator 4 Sd is 
applied continuously to the arithmetic unit. The clutter- 
lock section 46 is similarly controlled by the same tim- 
ing T. The arithmetic unit continuously applies these 
error signals Afo and Afb to arrive at the final parame- 10 
ters A, B, C and D as described above. The <f>(u) values 
are used to look up complex azimuth reference coeffici- 
ents for correlation, and the value <j>(t) divided by bi, 
and then divided again by b 2 , is used for the respective 
purposes of looking up interpolation weights and delay 15 
coefficients. 

All of these quantities are first stored in a reference 
storage unit 50 and are read out according to an FFT bit 
- reversal sequence required to supply the transversal 
filter with the appropriate coefficients according to the 20 
flow diagram of FIG. 10. Note that the 8 points (each of 
which is a line of di samples along the range direction 
for this exemplary diX8 data block are numbered in 
sequence on the left from 0 through 7, but after process- 
ing through the block-pipelined FFT, the order is not 25 
sequential, but rather as shown on the right. However, 
upon observing the binary numbers in parenthesis for 
the points 0-7 with the binary numbers of the output 
samples directly across, it is seen that by simply revers- 
ing the order of the bits for the input sample numbers, 30 
the proper binary number for the output samples is 
generated. For example, the output number for the 
input sample number 4 is 1 generated by reversing the 
sequence of the binary digits 100 to 001. 

Although particular embodiments of the invention 35 
have been described and illustrated herein, it is recog- 
nized that modifications and variations may readily 
occur to those skilled in the art. Consequently, it is 
intended that the claims be interpreted to cover such 
modifications and equivalents. 40 

What is claimed is: 

1. In a synthetic aperture radar system (SAR) having 
means for correlating range sampled data with a range 
reference function to produce the SAR response of 
range correlated data in sequential form, range line by 45 
range line in the time domain, a method for correlating 
said SAR response in the azimuth direction utilizing a 
block-pipelined Fast Fourier Transform (FFT) opera- 
tion with a predetermined FFT transform size com- 
prised of the step of delaying said SAR response to 50 
embed in the Fourier transform operation a corner-turn- 
ing function as said SAR response is converted from the 
time domain to a frequency domain, providing prede- 
termined delays for said SAR response for accurate 
range migration compensation and multiplying said 55 
SAR response in the frequency domain with an azimuth 
reference function in the azimuth frequency domain. 

2. A method as defined in claim 1 further including an 
inverse FFT operation to restore azimuth correlated 
data in the frequency domain to the time domain for 60 
imaging. 

3. A method as defined in claim 2 including the step 
of generating said azimuth reference function from 
estimates of Doppler center frequency, fa and fre- 
quency rate, f b, and including the step of continually 65 
feeding back single-look output data from said inverse 
FFT for adjustments in fo and fo estimates, whereby 
clutterlock and autofocusing is provided. 


357 

16 

4. A method as defined in claim 3 including the step 
of storing said azimuth reference function to be applied 
to sequential azimuth data points, where each point is a 
line of range samples, and reading out said azimuth 
reference function in a different order required to sup- 
ply said transversal filter with reference coefficients in 
appropriate order for block-pipelined FFT operation. 

5. A method as defined in claim 4 wherein said azi- 
muth data points are in binary numbered sequence, and 
said different order for reading out said reference coeffi- 
cients is established by reversing the sequence of binary 
digits for each number of data points in sequence and 
reading out the reference function stored for the binary 
number corresponding to the reverse binary digit for 
each number of data points in sequence. 

6. A method as defined in claim 5 including the step 
of generating for said transversal filter a range migra- 
tion compensation function, D, and interpolation 
weights, W i, as, a function of said frequency, f o, and 
frequency rate, fo, for delaying data into said transver- 
sal filter, and said interpolation weights being applied to 
i data points combined into a single output data point, 
where each data point is a line of range samples. 

7. A method as defined in claim 6 including the step 
of storing said delay function D and interpolation 
weights W / to be applied to sequential azimuth data 
points, where each point is a line of range samples, and 
the further step of reading out said azimuth reference 
function in a different order required to supply said 
transversal filter with said delay function D and inter- 
polation weights W ,• in the same bit-reversed sequence 
as said reference coefficients. 

8. In a synthetic aperture radar system (SAR) having 
means for correlating range sampled data with a range 
reference function to produce the SAR response of 
range correlated data in sequential form, range line by 
range line in the time domain, apparatus for correlating 
said SAR response in the azimuth direction utilizing a 
block-pipelined Fast Fourier Transform (FFT) having a 
predetermined FFT transform size comprised of means 
for delaying said SAR response to embed in the Fourier 
transform operation a corner-turning function as said 
SAR response is converted from the time domain to a 
frequency domain, and including a transversal filter 
connected to receive said SAR response in the fre- 
quency domain, said transversal filter having means for 
accurate range migration compensation and means for 
multiplying said SAR response in the frequency domain 
with an azimuth reference function in the azimuth fre- 
quency domain. 

9. Apparatus as defined in claim 8 further including a 
block-pipelined inverse FFT to restore azimuth corre- 
lated data in the frequency domain to the time domain 
for imaging. 

10. Apparatus as defined in claim 9 including means 
for generating said azimuth reference function from 
estimates of Doppler center frequency, fo, and fre- 
quency rate, fo, and including means for continually 
feeding back single-look output data from said inverse 
FFT to said generating means for adjustments in f/jand 
fo estimates, whereby clutterlock and autofocusing is 
provided. 

11. Apparatus as defined in claim 10 including means 
for storing said azimuth reference function to be applied 
to sequential azimuth data points, where each point is a 
line of range samples, and means for reading out said 
azimuth reference function in a different order required 
to supply said transversal filter with reference coeffici- 
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ents in appropriate order for block-pipelined FFT oper- 
ation. 

12. Apparatus as defined in claim 11 wherein said 
azimuth data points are in binary numbered sequence, 
and said different order for reading out said reference 5 
coefficients is established by means for reversing the 
sequence of binary digits for each number of data points 
in sequence and means for reading out the reference 
function stored for the binary number corresponding to 
the reverse binary digit for each number of data points 
in sequence. 

13. Apparatus as defined in claim 12 including means 
for generating a range migration compensation func- 
tion, D, as a function of said frequency, fa and fre- 
quency rate, fa and interpolation weights, W,-, for said 15 
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transverse filter, means responsive to said function D 
for delaying data into said transversal filter, means for 
applying said interpolation weights to i data points com- 
bined into a single data point, where each data point is 
a line of range samples. 

14. Apparatus as defined in claim 13 including means 
for storing said delay function D and interpolation 
weights W / to be applied to sequential azimuth data 
points, where each point is a line of range samples, and 
10 means for reading out said azimuth reference function in 
a different order required to supply said transversal 
filter with said delay function D and interpolation 
weights W,- in the same bit-reversed sequence as said 
reference coefficients. 
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